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ABSTRACT 

We present spectroscopic observations of ultra compact dwarf (UCD) galaxies in the 
Fornax and Virgo Clusters made to measure and compare their stellar populations. 
The spectra were obtained on the Gemini-North (Virgo) and Gemini-South (Fornax) 
Telescopes using the respective Gemini Multi-Object Spectrographs. 
■ We estimated the ages, metallicities and abundances of the objects from mea- 

surements of Lick line-strength indices in the spectra; we also estimated the ages and 
metallicities independently using a direct spectral fitting technique. Both methods re- 
vealed that the UCDs are old (mean age 10.8 ± 0.7 Gyr) and (generally) metal-rich 
(mean [Fc/H] = —0.8 ± 0.1). The alpha-element abundances of the objects measured 
from the Lick indices are super-Solar. 

We used these measurements to test the hypothesis that UCDs are formed by the 
tidal disruption of present-day nucleated dwarf elliptical galaxies. The data are not 
consistent with this hypothesis because both the ages and abundances are significantly 
higher than those of observed dwarf galaxy nuclei (this docs not exclude disruption of 
an earlier generation of dwarf galaxies) . They are more consistent with the properties 
of globular star clusters, although at higher mean mctallicity. The UCDs display a 
very wide range of metallicity (—1.7 <[Fc/H]< 0.0), spanning the full range of both 
globular clusters and dwarf galaxy nuclei. 

We confirm previous reports that most UCDs have high metalliticities for their 
luminosities, lying significantly above the canonical metallicitiy-luminosity relation 
followed by early-type galaxies. In contrast to previous work we find that there is no 
significant difference in either the mean ages or the mean metallicities of the Virgo 
and Fornax UCD populations. 

Key words: galaxies: formation - galaxies: star clusters: general - galaxies: stellar 
content - galaxies: dwarf 



1 INTRODUCTION 

Ultracompact dwarf (UCD) galaxies were originally discov- 
ered as very compact stellar sys tems (C SS) in the nearby 
Forn ax galaxy cluster (|Hilker et al. . 1999: Drinkwater et al.l . 
l200d ). Although as luminous as dwarf galaxies (—11 < 
Mb < —12), these objects were easily confused with Galac- 
tic stars in ground-based imaging. Measurements of the 
structural properties of these objects established that they 
were distinct from both globular clusters (more luminous 
and more distant from galaxies), and dwarf galaxies (much 
more compact) (|Phillipps et al. I, 1200 ll : iDrinkwater" et al.l . 
I2OO3I ). UCDs have since been found in the Virgo cluster 



E-mail: katy.francis@uqconnect.cdu.au 



( Jones et al. . 200a) and ot her more distant clusters (see 
iBaumgardt fc Mieskd . [2008I . for examples). 

The location of the UCD populations in dense galaxy 
clusters close to central galaxies with rich globular clus- 
ter populations prompted two main hypotheses for their 
formation. One model is that UCDs formed when the 
outer layers of nucleated dwarf galaxies are removed 
by t id al str ip ping near the cluster centre |Bekki et aT 



by t id ai str ipp 
200ll. I2OO3I; L 



iThomas. Drinkwater fc Evstigneeval . I200S 



Paudel. Lisker fc Janj . I2OI0I ). The other model is that 



UCDs are simply the extremely bright tail of the 
usual globular cluster (GC) pop ulations around t h e cen- 
tral galaxies in each c l uster (|Evstigneeva et aL 

I l2007l : 

iMieske. Milker fc Misgeldl . |2012| ). It is not as easy as first 
thought to separate these two models, so the question 
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of UCD fo rmation remains a topi c of active investi- 
gation (e.g. Chilingarian et al. . 2011 : Brodie et all I2OIII : 



IChiboucas et all I2OIII : iBriins et all . I2OIII ). 

Many studies have measured the structural properties 
of UCDs, benefiting in particular from Hubble Space Tele- 
scope imaging to resolve the objects. The latest analyses of 
the internal dynamics (as well as their stellar mass func- 
tions) show no need for dark matter in these objects to ex- 
plain their mass-to-light ratios (|Baumgardt fc Mieskd . [2OO8I : 
IChihngarian et all. [20T1I ') in contrast to some earlier results. 
On the other hand, there is increasing evidence that the 
more massive UCDs exhibit structural relationships unlike 
GCs. UCDs more luminous than Mv ~ —11 display a 
strong correlation between size and luminosity, whereas the 
sizes of globu l ar clu s ters are not correlated w ith luminosity 
iMieske et all (|2006l '): lEvstigneeva et all (|2008l '). 

The other major observational property of UCDs — 
the focus of this paper — is their chemical composition. 
The availability of 8-m class telescopes makes it quite 
feasible to make spectroscopic estimates of the ages and 
metallicities of UCDs in the nearer clusters. iMieske et al.l 
(j2006f l obtained Magellan Telescope spectra of 26 com- 
pact Fornax objects: the bright UCDs had high metallic- 
ities compared to dE galaxies in Fornax, suggesting they 
were very young, although the ages were not measured di- 
rectly. In c ontrast Keck spect r a of a sample of Virgo clus- 
ter UCDs jEvstigneeva et all . |2007| ) showed they were all 
old (more than 8 Gyr) like globular clusters. More evidence 
for a difference between the two cluster s was presented by 
iFirth. Evstigneeva fc Drinkwateil (|2009l ) who also reported 
young ages for some Fornax cluster UCDs, although these 
were measured wi t h low er quality spectra. More recently, 
IChilingarian et al.1 (|201ll ) analysed VLT spectra of 24 For- 
nax UCDs finding old ages consistent with globular clusters. 

In this paper we use high-quality data to make accurate 
measurements of the chemical composition of UCDs in both 
the Fornax and Virgo clusters. Our aim is to provide further 
tests of the formation models discussed above. In particular 
we wish to test the previous claims that the Fornax UCDs 
are younger than those in Virgo. We will do this by using 
data from similar instruments and identical analysis for both 
clusters to minimise any systematic effects. 

We describe our observational data and spectral analy- 
sis in Section[2] This consists of high signal-to-noise long slit 
spectra of Fornax and Virgo UCDs obtained with the Gem- 
ini North and South Muhi-Object Spectrographs (GMOS-N, 
and GMOS-S). In Section[3]we compare the results to stellar 
population models to determine the ages, metallicities, and 
a-element abundances of the UCDs. We discuss the impli- 
cations of the measurements in Section [4] and we summarise 
the paper in Section [S] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Selection of Targets 

We selected targets for our Gemini ob servations fr o m our 
larger UC D samples in the Fornax (ICregg et al.l. 20091) 
and Virgo (jJones et al. I I2OO6I : I Firth. Drinkwater fc Karickl . 
[2008. ) clusters. The targets were selected at random from 
the central (1 degree radius corresponding to a projected 



separation from the central galaxy of 200 kpc) field in 
each cluster, chosen so as to have luminosities in the range 
— 10 < Mb < —12. The only exception to this selection 
was that two Virgo objects from the "Intra-cluster field" 
were included in the Virgo sample. The mean luminosities 
of the targets were Mb = -11.0 ± 0.2 and -11.1 ± 0.2 for 
the Fornax and Virgo objects respectively^. There were no 
significant differences in the luminosity distributions in the 
two clusters. The UCDs measured are listed in Table[T]along 
with one globular cluster in each region observed for com- 
parison purposes. The table lists names for the objects based 
on their co ordinates except when n ame s were previously a s- 
signed by (|Drinkwater et al.l . [ioooh or (| Jones et al.U2006l ). 

For absolute magnitudes, we converted the SDSS g 
magnitudes to Cousins B magnitudes using the relation 
B = g + 0.42(g - r) - 0.024 (|Fukugita et al.l , i 19961 ). 

We adopt the following distance moduli in this paper: 
the Fornax Cluster at 31.35 mag, the Virgo Cluster at 30.97 
mag (see Table 1 of lFirth. Drinkwater fc Karickl . |200^ ) and 
the comparison globular cluster i n NGC 1407 at 31.99 mag 
jjerien. TuUv fc Trenthaml . l2004l ). 



2.2 Gemini Observations 

The spectroscopic data for ultra-compact dwarf galaxies 
were collected using GMOS-N and GMOS-S spectrographs 
on the 8-m GEMINI-North and GEMINI-South telescopes 
for Virgo and Fornax objects correspondingly. All obser- 
vations were obtained in the service mode by the GEM- 
INI observatory staff. Both spectrographs are equipped with 
6.3kx4.5k mosaic detectors comprising three 2kx4.5k CCD 
chips with narrow 90-pixel gaps between them. We de- 
scribe our processing of the two sets of data separately 
due to small differences in the setup. Data from one of 
the ob serving runs (GN-2007A-Q76) we r e pre viously pub- 
lished (|Firth. Evstigneeva fc Drinkwaten . l2009h , but we re- 
processed those data for the current paper to obtain a uni- 
form set of measurements. 



2.2.1 Gemini-North Virgo Data 

The data for Virgo cluster objects were collected for pro- 
posals GN-2007A-Q-76 (P.I.: Firth) and GN-2008A-Q-92 
(P.I.: Evstigneeva) in February-June 2007 and February- 
May 2008. We used the long-slit setup (slit length 5.5 ar- 
cmin) of GMOS-N with the B600-f _G5303 grating provid- 
ing a wavelength coverage between 3850 and 7050A. The 
slit contains two narrow bridges splitting it into three equal 
pieces used to trace and correct the geometric distortions. 
The binning of 2x4 yielded a spatial scale of 0.29 arc- 
sec pix~^ and a dispersion of «lA pix~^. The inter-chip gaps 
in our setup were around 4500A and 5760A. The slit width 
was set to 0.75 arcsec for VUCD 1232-1-0944 and 1 arcsec 
for all the remaining targets. Three 2300-sec long exposures 
were obtained for every object. Bias, flat and copper-argon 
arc line frames were obtained as a part of the standard cali- 
bration plan for GMOS-N. Beside scientific targets, spectra 
of two spectrophotometric standard stars, Hiltner 600 and 
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Table 1. Summary of the Targets 



Name(s) 


Run" 


Slit Width 




RA 


Dec 


Velocity 


Magnitude 


Radius^ 






^arcsec^ 




(h.m.s) 


\ ■ ■ } 


(km 


s-i) 


(mag) 




(arcmin) 


r UUUUooD-OOOD 


3 


0.75 


4800 


Uo;od:zz.zo 


-oo:od:o4.3 


1297 


± 22 


20.29 


r 

19.77 


27.50 


FUCD0336-3522 


3 


0.75 


4800 


03:36:26.72 


-35:22:01.6 


1290 


± 45 


20.20 


19.42 


25.41 


FUCD0336-3514 


3 


0.75 


4200 


03:36:27.74 


-35:14:13.9 


1416 


± 14 


20.12 


19.27 


27.84 


FUCD0336-3548 


3 


0.75 


7500 


03:36:47.74 


-35:48:34.1 


1416 


± 18 


20.74 


20.04 


29.80 


UCDl 


3 


0.75 


3000 


03:37:03.30 


-35:38:04.6 


1563 


± 14 


19.79 


19.03 


20.65 


FUCD0337-3536 


3 


0.75 


4800 


03:37:24.91 


-35:36:09.7 


1502 


± 35 


20.26 


19.61 


15.93 


FUCD0337-3515 


3 


0.75 


7500 


03:37:43.56 


-35:15:09.6 


1283 


± 21 


20.79 


20.12 


15.05 


FUCD0338-3513 


3 


0.75 


4200 


03:38:23.75 


-35:13:49.5 


1565 


± 25 


20.13 


19.57 


13.23 


FUCD0339-3519 


3 


0.75 


4200 


03:39:20.56 


-35:19:14.6 


1466 


± 13 


20.17 


19.42 


13.07 


UCD5 


3 


0.75 


3000 


03:39:52.58 


-35:04:24.1 


1282 


± 18 


19.86 


19.25 


28.33 


NGC1407-GC1 


3 


0.75 


7500 


03:40:09.42 


-18:33:37.3 


1219 


± 27 


B = 20.98= / 


= 19.16'= 


1.33 



VUCD1230-I-1233 


2 


1.0 


3100 


12:30:47.40 


VUCD2 


1 


1.0 


7500 


12:30:48.24 


VUCD3 


2 


1.0 


3900 


12:30:57.40 


Strom 417 


1 


1.0 


5100 


12:31:01.29 


VUCD1231-I-1234 


2 


1.0 


4800 


12:31:02.59 


VUCD5 


2 


1.0 


5400 


12:31:11.90 


VUCD8 


1 


1.0 


7200 


12:32:04.33 


VUCD9 


1 


1.0 


6900 


12:32:14.61 


VUCD1232-I-0944 


1 


0.75 


6900 


12:32:55.55 


VUCD 1233-1-0952 


2 


1.0 


5100 


12:33:07.36 











g 


/ 


r / 




12:33:01.70 


1390 


± 


20 


19, 


.53 


18.89 


9.57 


12:35:11.10 


886 


± 


19 


19, 


.21 


18.54 


11.72 


12:25:44.80 


640 


± 


11 


18, 


.87 


18.07 


3.00 


12:19:25.60 


1872 


± 


15 


19, 


.69 


19.01 


4.97 


12:34:14.11 


1198 


± 


23 


19, 


,81 


19.21 


11.24 


12:41:01.20 


1287 


± 


13 


19, 


,10 


18.39 


18.39 


12:20:30.62 


1679 


± 


18 


19, 


.63 


19.05 


18.54 


12:03:05.40 


1346 


± 


16 


19, 


.41 


18.84 


29.13 


09:44:12.70 


1699 


± 


37 


20, 


.84 


20.26 


111.14 


09:52:54.30 


1114 


± 


18 


19, 


.97 


19.39 


123.27 



"Gemini project codes: 1 = GN-2007A-Q76, 2 = GN-2008A-Q-92, 3 = GS-2008B-Q-6 
''Total exposure time 

'^Projected radius from nearest large galaxy 
''CCD Photometry (Karick in preparation) 

Photometry by Cenarro et al. (2007) 
•''SDSS photometry fFirth. Drinkwater fc Karick. 20081 



HZ 44 were obtained in order to perform the flux calibra- 
tion, as well as the twilight spectral frames which we used 
to estimate the spectral line spread of GMOS-N. 

We reduced the data using our own GMOS data reduc- 
tion pipeline constructed on top of the universal IFU data 
reduction toolbox implemented in idl. The data reduction 
was done independently for every science exposure (includ- 
ing standard stars and twilight flats). 

The data reduction involved several stages. The flrst 
was executed for every individual chip and included the fol- 
lowing steps: 

• Bias subtraction, masking bad columns and hot pixels, 
removi ng cosmic ray hits u sing the Laplacian filtering tech- 
nique (|van Dokkumll200ll '). 

• Correcting counts for the read-out analog-to-digital 
gain. 

• Modelling the global diffuse light by low-order two- 
dimensional polynomial surface using counts outside the slit 
(about 15 binned pixels below and above it) and slit bridges 
and subtracting it. 

Then, in order to mosaic of individual chips, we used the 
transformation coefficients for relative offsets and rotations 
determined from the IFU arc line and fiat frames obtained 
for a different observing programme and found in the GEM- 



INI science archivfl IFU data frames obtained without the 
detector binning contain traces of individual fibers. They al- 
lowed us to determine the transformation coefficients with 
much higher precision that can be achieved using our binned 
frames in the long-slit mode. 

We performed the following data reduction steps on the 
mosaiced frames: 

• The normalised fiat field was created by removing the 
overall continuum shape by dividing every line by the mean 
flat field vector along the slit; then it was used to correct all 
science and calibration frames. 

• Arc lines were identified using the line list from the 
official IRAF GMOS reduction pipeline and an iterative pro- 
cedure was usedto remove blended lines and those heavily 
affected by inter-chip gaps and bad columns on individual 
detectors. Then, the two-dimensional dispersion relation was 
built using the two-dimensional polynomial surface of the 
5rd and 3rd orders along and across the wavelength range 
respectively. The RMS of the residuals of the dispersion re- 
lation was usually about 0.09 pix. 

• We created a model of the ni ght sky airglow spectrum 
using an approach proposed by iKelsonI (|2003l ) aimed at 
avoiding the artefacts created by the interpolation of nar- 
row airglow lines. We used the non-linearised science frame 
and a two-dimensional dispersion relation in order to con- 

^ http:/ /www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/gsa/ 
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struct an oversampled spectrum of the airglow emission in 
the wavelength space by using the slit regions not contami- 
nated by the scientific targets or other serendipitous targets. 
This spectrum was then approximated using a smoothing b- 
spline parametrization with equidistant nodes every O.SA. 
Then, this parametrization was used to evaluate the airglow 
spectrum at every position on the slit again using the two- 
dimensional wavelength relation. The sky subtraction made 
in this fashion has nearly Poisson quality. 

• After the sky subtraction, we linearised the spectra. 

• Then, we performed the flux calibration using the spec- 
tra of spectrophotometric standard stars. 

The flnal step was the optimal extraction of the source 
on the slit which is made using an empirically constructed 
point-spread function from the spectrum itself as all our tar- 
gets are sufficiently bright. Then, the spectra from individual 
exposures were combined. Since for every target the individ- 
ual exposures were obtained during rather short periods of 
time (an order of a few days), the individual heliocentric 
corrections were not required. 

The uncertainty frames were computed from the pho- 
ton statistics and the read-out noise values, and processed 
through exactly the same data reduction steps in order to 
estimate the flux uncertainties. 

For certain targets, arc line frames were not obtained 
for every science exposure, however by comparing the indi- 
vidual science frames, notable shifts (up-to 2-3 pixels) have 
been detected both along and across dispersion. For such 
situations, we took one of the available arc line frames and 
shifted it according to the measured offsets. We stress, that 
one should not apply shifts to science exposures as they 
would introduce additional artefacts because of the interpo- 
lation. 



2.2.2 Gemini-South Fornax Data 

The data for Fornax cluster objects were collected for 
proposal GS-2008B-Q-6 (P.I.: Drinkwater) in August- 
November 2008. As for the Virgo cluster galaxies, we used 
the long-slit setup (slit length 5.5 arcmin) of GMOS-S with 
a similar B600-|-_G5323 grating providing a wavelength cov- 
erage between 3900 and 7000A. As in GMOS-N, the slit con- 
tains two narrow bridges splitting it into three equal pieces 
used to trace and correct the geometric distortions. Simi- 
larly, we used binning of 2x4 that yielded a spatial scale of 
0.29 arcsec pix^^ and a dispersion of «lA pix^^. The slit 
width was set to 0.75 arcsec for all targets. An important 
difference with our GMOS-N observations was the spectral 
dithering used to flU the inter-chip gaps: the central wave- 
length were different by 50A for every of the three individual 
exposures per target. The GMOS-S calibrations were simi- 
lar to those obtained with GMOS-N, only the spectroscopic 
standard star was different (EG 131). 

We reduced the data using the same software as for 
GMOS-N data with slight modiflcations in the two steps. 

The arc line spectra obtained during the night time for 
corresponding observing blocks were read out only in the 
central part of the slit of about 1 arcmin in size (trun- 
cated arc line frames hereafter). This setup was sufficient 
to achieve high-quality wavelength calibration for the sci- 
ence spectra because all our objects were point sources cen- 



tered on the slit. However, the algorithm that we use for the 
sky subtraction required much longer slit as it exploited the 
curvature of spectral lines in the non-linearised CCD frames. 
Therefore, we modelled the missing regions of arc line frames 
for every observing block using full-frame arc line spectra 
obtained at the end of the observing programme. 

For our modelling we assumed that between individ- 
ual observing blocks: (1) slight overall shifts of the spectra 
of an order of a few pixels in both directions may present; 
(2) slight differences between the wavelength solutions may 
present and they could be precisely accounted by the 2nd or- 
der polynomial; (3) there was no significant rotation; (4) the 
distortion along the slit remained stable. The relative overall 
shift across dispersion was determined by cross-correlating 
the spectral flat fleld frames. The relative position of the slit 
bridges proved our 4th assumption. The shift along the dis- 
persion was accounted together with the differences between 
the wavelength solutions by individually cross correlating 
three fragments of the truncated and full arc line frames, 
near the centre and at the ends of the spectral range and 
then computing the 2nd degree polynomial coefficients de- 
scribing the differences. Then, the full-frame arc line spectra 
were transformed in order to match these shift using the 2D 
polynomial image warping and used to fill the missing data 
in the truncated arc line frames. 

Subsequent tests showed that assumptions 2 (that a 2nd 
order polynomial was sufficient) and 3 (that there was no 
significant rotation) were justified. Firstly, we did not de- 
tect any measurable shifts between the truncated arc line 
frames and the corresponding parts of the modelled frames. 
Secondly, we did not have any sky subtraction artefacts that 
would be obvious in case of systematic errors in the wave- 
length calibration derived from the modelled frames. 

Another difference with the GMOS-N reduction arose 
from the slightly different observing strategy: the spectral 
dithering with a pattern of 50A was applied in order to fill 
the inter-chip gaps. However, the spectrophotometric stan- 
dard star was observed only in one position. Therefore, due 
to the slight change of the grism blaze function as a result 
of the central wavelength shift, and because of the slight 
differences between the sensitivity curves of individual CCD 
chips constituting the GMOS-S mosaic detector, the quan- 
tum efficiency curve obtained from the spectrum of EG 131 
cannot be directly applied to calibrate the spectra obtained 
with different grism positions. We stress that these effects 
cannot be accounted by fiat fielding, because the continuum 
shape of the fiat field lamp was not flat. 

In order to perform the flux calibration of the spectra 
obtained using grism positions different from that of the 
EG 131 observation, we considered that the quantum effi- 
ciency curve could be represented as a low-order (9th) poly- 
nomial function with two breaks corresponding to the chip 
gaps, i.e. adding two multiplicative coefficients for the 2nd 
and 3rd CCD chips. We performed a non-linear minimiza- 
tion in order to determine these coefficients while the poly- 
nomial coefficients were fitted linearly at every evaluation 
of the function. The obtained values suggested the differ- 
ences of an order of 5-8 per cent between individual chips 
that we took into account. This approach was proved to be 
valid when we co-added individual spectra for UCDs since 
no breaks were detected (either visually or in the residuals 
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of the full spectral fitting) in the spectral regions around the 
CCD gaps of individual observing blocks. 

2.3 Measurement of Spectra 

We used two approaches here to estimate the ages and chem- 
ical compositions of the UCDs from the Gemini spectra. 
Following much of the previous analysis we started by mea- 
suring individual (Lick) spectral line indices, for comparison 
to stellar population models. We also applied a more gen- 
eral approach which derived the parameters by direct fit- 
ting of population models to the whole spectra. Due to the 
relatively low spectral resolution in the selected setups of 
GMOS-N and GMOS-S, velocity dispersion measurements 
were not possible. 

2.3.1 Lick Index Analysis 

In order to measure the Lick indices it was necessary to 
smoot h our spectra to the Lic k resolution of (JLick ~ 3.6 
A, fsee lEvstigneeva et akllioOTl V This was achieved by con- 
volving the spectra with Gaussians (cr = 3.03 A for the 1 
arcsec slit, and a = 3.26 A for the 0.75 arcs ec slit). We 
then used the indexf prograrrffl jCardie 1 120071') to measure 
the Li ck/IDS line strength indices, defined bv lWorthev et all 
(|l994 ). of each object. The uncertainties in the measure- 
ments were estimated using the uncertainty arrays described 
above, scaled appropriately for the smoothed spectra. The 
measured indices, along with the two derived indices (Fe) 
and [MgFe]' are listed in Table [2] The two derived indices 
are found as follows: (Fe) — average of Fe5270 and Fe5335, 
and [MgFe]' = [Mg&(0.72xFe5270 + 0.28xFe5335)]^''^ 
jThomas. Maraston fc Bended . [ioO^ I . 

The Lick indices offer the advantage of a substantial 
body of previous measurements for comparison and discus- 
sion, but they do not make use of all the information avail- 
able in the spectral data. One way to improve the analysis 
is to combine the values of mul tiple Lick indices to obtain 
more accurate results, as done bv lProctor. Forbes fc BeaslevI 
(|2004l ). In Section r2.3.2l we take a slightly different approach, 
by fitting models directly to the full spectra themselves 
rather than first calculating line indices. 



2.3.2 Full spectral fitting 



We used the nb ursts 
Chilingarian et al.l, 



full spectral 

2007bl lal) with 



fitting technique 
the pegase.hr 



Le Borgne et al.l . |2004| ) simple stellar population (SSP) 



models in order to determine SSP-equivalent ages and 
metallicities. 

The NBURSTS full spectral fitting package implements a 
pixel-space fitting algorithm. Generally, an observed spec- 
trum is approximated by a linear combination of stellar 
population models broadened with the galaxy's paramet- 
ric line-of-sight velocity distribution, whose parameters (e.g. 
age, metallicity, initial mass function) are determined in- 
side the same minimization loop as the internal kinemat- 
ics. The fitting procedure includes a multiplicative polyno- 
mial continuum aimed at absorbing possible flux calibra- 

^ |http : // www, ucm . es / info / Astrof / software / indexf / indexf. ht ml | 
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tion issues both, in observations and in the models. The 
NBURSTS technique is shown to produce unbiased esimates 
of ages and metallicities f or a-enhanced stellar populations 
I Chilingarian et all . l2008f ) although the pegase.hr steUar 
population models are not representative of such popu- 
lations because they are constructed from empirical stel- 
lar spectra for stars mostly from the Solar neighbourhood 
which are known to have their a/Fe abundance ratios anti- 
correlated with metallicities. 

In this study we use the simplest case of the fitting 
procedure, when the observed spectrum is approximated 
by one SSP model characterised by its age and metallicity. 
We used an SSP grid computed with the pegase.hr evolu- 
tionary synthesis code at the intermediate spectral resolu- 
tion (R=10000) in a wave l eneth range 3900-6800 A for the 
iKroupaT Tout fc Giliiior^ (| 19931 ) stellar initial mass func- 
tion. We used the 29th order multiplicative polynomial con- 
tinuum. The model SSP grid was pre-convolved with spec- 
tral line spread functions of GMOS-N and GMOS-S spectro- 
graphs varying with the wavelength. We computed them by 
fitting the twilight spectra obtained in the same instrumen- 
tal setups as UCD spectra against the R=10000 Solar spec- 
trum from the elodie.3 (|Prugniel fc Soubiranl . |2004 ) stel- 
lar library in five wavelength segm ents between 3900 < A < 
6800 A using the ppxf procedure (|Cappellari fc Emselleml . 
|2004 ). For the 1-arcsec wide slit of GMOS-N, the spectral 
line spread changes from crinst ~ 150 km s~^ at 3900 A to 
100 km s~^ at 6700 A, w ith nearly constant Gauss-Herm ite 
coefficients /i3 and /i4 (Ivan der Marel fc Franxl. 1 19931 ) of 
about 0.00 and —0.10 correspondingly. For the 0.75-arcsec 
wide slit of GMOS-S, the spectral resolution was slightly 
better and corresponded to ainst = 120 km s^^ at 3900 A 
and 70 km s"^ at 6700 A with the h3 and /i4 values of 0.00 
and -0.08. 

In addition to the best-fitting values of radial ve- 
locity and SSP-equivalent age and metallicity (included 
in Table ^ which are obtained in a single mini- 
mization loop, for every object we computed a con- 
fidence map in the age-metallicity space simil ar to 

those presented in JChilingarian. Cavatte fc Bergondl . l2008l : 
IChilingarian et all. [2OIII ). These can be found in Ap- 
pendix [X] We notice, that when the NBuRSTS technique 
is used in the x^-ma-pping mode, i.e. fitting a single SSP 
with fixed age and metallicity at every grid node, the al- 
gorith m "degenerates" into the form totally equivalent to 
PPXF ICappellari fc Emselleml . 12004) ') with a single template 
spectrum. 



3 ANALYSIS 

In this section we measure the internal composition of the 
UCDs from the Gemini spectra. We first compare the values 
of individual line indices to models before directly fitting 
stellar population models to the spectra. 



3.1 a-element Abundance Ratios from Lick 
Indices 

The star formation history of stellar systems, such as dwarf 
galaxy nuclei and GCs, can be probed with Q-element abun- 
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Table 2. Lick/IDS Atomic and Derived Indices 



Target Lick/IDS atomic indices Derived Indices 

H/3 H5A H7A Mg h Fe5270 Fe5335 (Fe) [MgFe]' 

(A) (A) (A) (A) (A) (A) (A) (A) 



FUCD0336-3536 


2 


61 


± 





05 


3.3 


± 


0.1 


0.91 


± 


0. 


09 





73 


± 





05 


1 


12 


± 





05 





59 


± 





06 





86 


± 





04 





84 


± 


0.03 


FUCD0336-3522 


1 


57 


± 





05 


-0.4 


± 


0.1 


-2.0 


± 





1 


1 


83 


± 





05 


2 


05 


± 





06 


1 


07 


± 





07 


1 


56 


± 





05 


1 


80 


± 


0.03 


FUCD0336-3514 


1 


40 


± 





05 


-2.3 


± 


0.14 


-6.4 


± 





1 


4 


62 


± 





04 


2 


62 


± 





05 


2 


17 







06 


2 


40 


± 





04 


3 


39 


± 


0.03 


FUCD0336-3548 


1 


99 


± 





05 


1.1 


± 


0.1 


-2.4 


± 





1 


2 


31 


± 





05 


1 


83 


± 





06 


1 


45 


± 





06 


1 


64 


± 





04 


2 


00 


± 


0.03 


UCDl 


1 


68 


± 





04 


-1.3 


± 


0.1 


-5.31 


± 





09 


3 


81 


± 





04 


2 


15 


± 





05 


2 


13 


± 





05 


2 


14 


± 





04 


2 


86 


± 


0.03 


FUCD0337-3536 


2 


29 


± 





05 


-6.6 


± 


0.2 


-2.4 


± 





1 


1 


87 







05 


1 


92 


± 





05 


1 


07 


± 





06 


1 


50 


± 





04 


1 


77 


± 


0.03 


FUCD0337-3515 


1 


76 







05 


-0.1 




0.1 


-3.5 







1 


3 


50 


± 





05 


2 


45 


± 





05 


1 


79 


± 





06 


2 


12 


± 





04 


2 


82 


± 


0.03 


FUCD0338-3513 


2 


35 


± 





04 


1.4 


± 


0.1 


-0.96 


± 





08 


1 


24 


± 





05 


1 


18 


± 





05 


1 


09 


± 





06 


1 


14 


± 





04 


1 


20 


± 


0.03 


FUCD0339-3519 


1 


78 


± 





05 


-1.5 


± 


0.1 


-5.3 


± 





1 


3 


52 


± 





05 


1 


97 


± 





05 


1 


77 


± 





05 


1 


87 


± 





04 


2 


60 


± 


0.03 


UCD5 


2 


25 


± 





05 


1.4 


± 


0.1 


-1.5 


± 





1 


1 


69 


± 





05 


1 


40 


± 





06 


1 


24 


± 





06 


1 


32 


± 





04 


1 


51 


± 


0.03 


NGC1407GC1 


2 


15 


± 





06 


0.9 


± 


0.2 


-1.7 


± 





2 


2 


28 


± 





06 


1 


95 


± 





06 


1 


60 


± 





07 


1 


78 


± 





05 


2 


05 


± 


0.04 


VUCD1230+1233 


1 


89 


± 





05 


-1.6 


± 


0.1 


-2.8 


± 





1 


2 


93 


± 





05 


2 


12 


± 





05 


1 


81 


± 





05 


1 


97 


± 





04 


2 


44 


± 


0.03 


VUCD2 


2 


21 


± 





04 


0.0 


± 


0.1 


-2.0 


± 





1 


1 


72 


± 





04 


1 


69 


± 





04 


1 


28 


± 





04 


1 


49 


± 





03 


1 


65 


± 


0.03 


VUCD3 


1 


40 


± 





04 


-5.6 


± 


0.1 


-7.6 


± 





1 


5 


10 


± 





03 


2 


55 


± 





04 


2 


22 







04 


2 


39 


± 





03 


3 


54 


± 


0.02 


Strom 417 


2 


04 


± 





04 


-0.9 


± 


0.1 


-1.2 


± 





1 


2 


82 


± 





04 


2 


31 


± 





04 


1 


88 


± 





04 


2 


10 


± 





03 


2 


48 


± 


0.02 


VUCD1231+1234 


1 


98 


± 





06 


0.8 


± 


0.1 


-1.2 


± 





1 


1 


29 


± 





06 


1 


23 


± 





06 





51 


± 





07 





87 


± 





05 


1 


15 


± 


0.04 


VUCD5 


1 


73 


± 





03 


-1.69 


± 


0.09 


-4.71 


± 





07 


3 


34 


± 





03 


2 


21 


± 





03 


1 


81 


± 





04 


2 


01 


± 





03 


1 


15 


± 


0.04 


VUCD8 


2 


13 


± 





04 


0.6 


± 


0.1 


-2.22 


± 





09 


1 


47 


± 





04 


1 


65 


± 





04 


1 


10 


± 





05 


1 


38 


± 





03 


1 


48 


± 


0.03 


VUCD9 


2 


43 


± 





04 


-0.17 


± 


0.09 


-2.61 


± 





08 


2 


07 


± 





03 


1 


77 


± 





04 


1 


60 


± 





04 


1 


69 


± 





03 


1 


89 


± 


0.02 


VUCD1232+0944 


1.3 


± 





1 


-6.5 


± 


0.5 


-4.3 


± 





3 


1 


27 


± 





09 


2 


01 


± 





08 


1.5 


± 





1 


1 


76 


± 





06 


1 


54 


± 


0.06 


VUCD1233+0952 


2 


12 


± 





07 


-2.0 


± 


0.2 


-2.6 


± 





2 


2 


02 


± 





07 


2 


16 


± 





07 


1 


19 


± 





08 


1 


68 


± 





05 


1 


95 


± 


0.04 



dance ratios (e.g. Geha. Guhathakurta fc van der Marell . 
I2OO3I : iBrodie fc StradeJ . |2006| . respectively'). 

In simple closed models of stellar evoluti on (e.g. 
lMcWilliaml . ll997l : lBressan. Chiosi fc Fagotto! . 1 1994h . the av- 
erage metallicity ([Fe/H]) of the populations always in- 
creases with time as metals are formed in stars and ejected 
into the ISM by supernovae. The high-mass supernovae 
(type II) produce large amounts of the "a-elements" (such 
as oxygen and magnesium) by neutron capture processes 
leading to a high a-element abundance ([a/Fe]). After a de- 
lay (due to their lifetimes) the low-mass stars stars will start 
producing type la supernovae: these mainly produce iron, so 
the alpha abundance then starts to decrease. 

In a given observed stellar population, low (Solar 
or sub-Solar) a-abundances reflect star formation which 
has taken place at a low rate over a long period, so 
that the type la supernovae have significantly reduced 
the [a/Fe] r atio. The cores of dwarf galaxies sho w these 
properties 



( Geha. Guhathakurta . 
1997^. 



van der Maxell . l2003l : 



iGorgas et al.l . Il997f ). In contrast, high (super-Solar) a- 
element abundances reflect a rapid early stage of star 
formation which formed most of the population. This 
produces much larger amounts of the a?-elements be- 
fore the SN la start, so the [a/Fe] ratio remains high. 
Globular clusters typically have super-Solar abundances 
ICohen. Blakeslee fc Rvzhovl Il998l ) having formed all their 
stars very early in the Universe. 

We present measurements of the Lick indices sensi- 
tive to a-element abundance for the UCDs in Figure [T] 
The figure plots the derived iron index (Fe) (see pre- 
vious section) against Mgb. The figure also shows pre- 
dicted model grids for populations with a range of abun- 
a/Fe]= -0.3,0,0.3,0.5). The in odcls wore calcu- 
Thomas. Maraston fc JohanssonI |2011. ) . These new 



dance ( 
lated by 

models have been calculated for Lick indices measured 



with flux-calibrated spectra, so do not require any zero- 
point corrections (such as used i n our earlier analysis — 
iFirth. Evstigneeva fc Drinkwatetl . l2009h . 

The figure shows that the majority of the UCDs in both 
clusters have super-Solar a-element abundances, mostly 
consistent with the [a/Fe]= 0.3 model. However, some of the 
UCDs with higher (Fe) values are also consistent with the 
[a/Fe]= 0.5 model. These high abundances are consistent 
with both old elliptical galaxies (shown as crosses on the fig- 
ure) and glob ular clusters (open triangle and square, bu t also 
see Fig. 2 of [rii omas. Maraston fc JohanssonI . I2OIII ) . The 
average abundance of the UCD population (parameterised 
by the ratio (Fe) /Mgb) is significantly higher (p = 0.003 us- 
ing a T-test) than that of nucleated dwarf galaxies (shown 
as open pentagons in the figure). There is no significant dif- 
ference between the abundances of the Virgo and Fornax 
UCDs. We note that one UCD (VUCD1232-f 0944) has a 
lower abundance, consistent with the Solar abundances of 
the dE,N galaxies. 



3.2 Metallicity and Age from Lick Indices 

Several of the Lick indices are sensitive to metallicity and 
age. We chose to focus on the combination of II/3 and the 
compound ind ex [MgFe]' for ease of comparison with p revi- 
ous work (e.g. iFirth. Evstigneeva fc Drinkwaterl . |2009| V In 
addition, the [MgFe]' index has the advantage that it mea- 
sures metallicity without being affected by a-element abun- 
dance (|Thomas. Maraston fc Bendeil . lioOSh . 

We show these indices measured for the UCDs in Fig- 
ure [2I along with the predicted mode l grid for a SSP 
model (|Thomas. Maraston fc Johanssoia . l201lf ) with the 
super-Solar abundance ([a/Fe] ~ 0.3) found in the pre- 
vious section for these objects. The figure also shows 
for comparison purposes the positions of samples of gi- 
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[a/Fe] = -0.3 



0.0 



4 - 



fO.3 



2 - 




Fornax UCD, GC * A 

Virgo UCD. GC ■ □ 

pub Virgo UCD X 

pub. dEN, E O • 



Mgb (A) 



2 - 



Fornax UCD, GC A A 

Virgo UCD, GC ■ □ 

pub. Virgo UCD X 

pub. dEN, E O • 

pub. GC o 




Age (Gyr) 



2 3 
[MgFe]' (A) 



Figure 1. a-element abundances of UCDs based on Lick in- 
dex measurements. SSP model grids are plotted fo r a ra nge of 
[a/Fe] values from lThomas. Maraston &: JohanssonI l|201lf l. The 
grid lines simulate ranges of constant age and metallicity. The 
Gemini data are plotted as (red) triangles for the Fornax ob- 
jects and (blue) squares for the Virgo objects. In each case the 
solid symbols are UCDs and the ope n symbol is a GC. Pre - 
vious measurements of Virgo UCDs l lEvstigneeva et al.l. |2007^ 
are plotted as (magenta) crosses. Previously published data are 
shown for comparison: nucleated dE galaxies (unfilled pentagons 
ICeha. Guhathakurta fc van der Marei|2003l l and elliptical galax- 
ies (filled pentagons iTrager et all I2OO0I) ! 



ant elliptical galaxies, nucleated dwarf elliptical galax- 
ies, and globular clusters. The dE,N galaxies display 
the youngest ages, with slightly sub-Solar metallicities. 
The giant ellipticals are older, with higher metallici- 
ties, whereas the glob ular clusters are also old, bu t with 
low metallicities (see Cohen. Blakeslee fc Rvzhovl 1 19981 : 



iGeha. Guhathakurta fc van der Mareil2003l ). We note that 
the globular clusters in Figure [2] show a considerable scat- 
ter below the lowest (maximum age) of the model grid lines. 
This is a known effect, possibly related to the modellin g 
of hot giant stars (|Thomas. Maraston fc JohanssonI . I2OIII ). 
The higher order H7 index shows less scatter, as shown in 
Figure |3] (without the comparison objects), but the overall 
conclusions are not changed (i.e. the UCD positions com- 
pared to the models). 

We inverted the SSP models 

(|Thomas. Maraston fc JohanssonI . I2OI j ) plotted in Fig- 
ure [2] to obtain estimates of the ages and metallicities 
of ea ch object using a bin o mial polynomial int e rpola- 
tion (ISanchez-Blazguez et al.l . l2006al : ICardiel et all . l2003l : 
IWolberd . Il990t ). The corresponding ages and metallicities 
obtained from the Lick indices are given in Table [S] 

Both of the figures show the UCDs clustering along the 
oldest model track, but with a considerable range of metal- 
licity. There is a scatter below the oldest model track in 
the H/3 plot but this is consistent with the scatter observed 
for the globular clusters. This scatter in age leads to some 



Figure 2. Age and metallicity of UCDs as indicated by the H/3 
and [MgFe]' L ick indices. SSP model grids are p l otted for [o/Fe] = 
0.3 dex from |Th omas. Maraston &: JohanssonI l|20l3)- The grid 
lines arc simulated constant age and metallicity. The symbols are 
as in the previous figures with the addition of M87 (Vi r go Cl uster) 
globular clusters (circles [Cohen. Blakeslee &i Rvzhovl . Il998l) . 
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Figure 3. Age and metallicity of UCDs as indicated by the H7 
and [MgFe]' L ick indices. SSP model grids are p l otted for [o/Fe] = 
0.3 dex fromlrh omas. Maraston fc Johansso 3 1 120111) . The grid 
lines are simulated constant age and metallicity. The comparison 
objects from the previous plots are not shown. This plot shows 
that, although the H7 index does not display as much scatter as 
H/3, the UCD have similar properties compared to the SSP model. 
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overlap with the dE galaxies, but the mean age of the UCD 
sample (as parameterised by [MgFe]'/2+H/3) is significantly 
(p < 0.001 with a T-test) older than that of the dE popu- 
lation. As above, we find no evidence for any significant age 
difference between the Virgo and Fornax UCDs. The UCDs 
have a much wider range of metallicity than any of the com- 
parison galaxy populations: the range is however consistent 
with that of the c omparison GC population (from M 87 in 
the Virgo Cluster; ICohen. Blakeslee fc Rvzhovlll998l ). The 
UCD population therefore includes a metal-rich component 
like many GC populations. 

We also show so me previous UCD measu rements using 
the Keck Telescope (|Evstigneeva et all 120071 ') in Figures [T] 
and[2l The object with the highest metallicity from the Keck 
measurements is the same object we have repeated in our 
Gemini observations, VUCD3. In each case both measure- 
ments are consistent, giving an indication of the repeatabil- 
ity of our measurements. 

One object, VUCD1232-f 0944 in the Virgo Cluster, has 
extreme values in all the index plots. It has the lowest abun- 
dance of the UCDs (Solar) in Figure [T] and also the lowest 
H/3 index in Figure (2] The H-gamma value also shows this 
object lying below the main population. We inspected the 
spectrum of this object and did not find any unusual fea- 
tures. The population model gives an age of 12.8 Gyr and 
[Fe/H] = —1.0, which is consistent with the Lick results 
given the large scatter in the II/3 index. 

3.3 Metallicity and age from the full spectral 
fitting 

The NBuRSTS spectral fitting approach makes full use of 
the available wavelength range and spectral resolution of 
the data so it should provide more reliable ages and metal- 
licities than those derived from the line indices which only 
use a small fraction of the available age- and metallicity- 
sensitive data contained in the index definition regions. 
However, a comparison of the metallicity estimates in Fig- 
ure U reveals a systematic metallicity-dependent offset: the 
NBuRSTS metallicities are about 0.2 dex lower than the Lick 
values at high metallicity. 

This offset may be caused by the models used by 
NBuRSTS: the age estimates from the full spectral fitting are 
not affected by n on-solar abundance ratio s of stellar pop- 
ulations (see e.g. IChilingarian et all . |2008| ) but the metal- 
licities can be affected. This is a result of the properties of 
the pegase.hr stellar population models used by NBuRSTS. 
These models are based on an empirical stellar library con- 
taining only stars in the Solar vicinity (see above) hav- 
ing a characteristic abundance pattern, representative of a- 
enhanced ([Mg/Fe] ~ -fO.3 dex) populations at low iron 
metallicities ([Fe/H] < —1.0 dex) but Solar abundances 
([Mg/Fe] « 0.0 dex) at high, i.e. Solar and slightly sub- 
Solar, iron abundances. Most of the objects in our study 
are «-enhanced (see above). Although the spectral region 
contains many features of a-elements between 4200A and 
5200A, the fitting residuals of metal-rich UCDs such as 
VUCD3 exhibit a strong template mismatch at the Mg6 
triplet position. Therefore the fitting procedure is represen- 
tative of iron metallicities when fitting spectra a-enhanced 
metal-rich objects where models are Solar-scaled in [a/Fe]. 
At the same time, for metal-poor targets we are measuring 




[Z/H] from Lick indices (dex) 

Figure 4. Comparison of metallicities ealculated from the Lick 
indices and from the full spectral fitting. The theoretical line cor- 
responding to full agreement between the methods is indicated in 
black: there is an offset of about 0.2 dex as discussed in the text. 

the mean because the models themselves are representative 
of a-enhanced populations. This explains why at high metal- 
licities our [Z/H] are slightly underestimated (about 0.2 dex) 
compared to the abundance-insensitive measurements done 
with the Lick indices. 

Generally, ages are determined with much higher uncer- 
tainties (in dex) than for the metallicities. This illustrates 
that most features in the spectrum are much less sensitive to 
age than to metallicity. For this reason we have not made a 
detailed comparison of the ages derived by the two methods. 

Confidence level maps of stellar population parameters 
in the age-met allicity space demonstrate quite well a well- 
known age-metallicity degeneracy which is manifested by 
the inclined shape of error contours so that older ages may 
"compensate" lower metallicities. The degeneracy is much 
weaker though for metal-poor objects. Generally, ages are 
determined with a much higher uncertainties (in dex) than 
metallicities that illustrates the fact that most features in 
the spectrum are much less sensitive to the age than to the 
metallicity. 

In many cases (5 in Fornax and 5 in Virgo) we ob- 
serve double minima in the age-met allicity confidence level 
maps with one of them always at very old age values signifi- 
cantly exceeding the presently accepted age of the Universe 
(13.7 Gyr) and reaching the maximal age in the model grid 
(17 Gyr). Their origin is unclear and it is probably connected 
to the model imperfections in the selected wavelength range. 
Although in a few cases those "old" minima are primary (i.e. 

has a global minimum there) , we discard all age estimates 
exceeding 15 Gyr and consider the secondary minimum if it 
exists. 

Four of the objects from our present Fornax 
Cluster sample have been previously analysed us- 
ing the NBuRSTS technique applied to high-resolution 
spectra: FUCD0337-3538=UCD1, FUCD0339-3504=UCD5, 
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Table 3. Age and Motallicity results from Liek indiees and Full Spectral Fitting 



from Lick Indices: from Spectral Fitting: " 

Target Age Metallicity Age Metallicity 

(Myr) (dex) (Myr) (dex) 



FUCD 0336-3536 


7 980 


± 


200 


-1.58 


± 


0.04 


> 15 000 '' 


-1.76 


± 


0.02 


FUCD0336-3522 


11 030 


± 


160 


-0.70 


± 


0.01 


8 180 ± 230 


-0.82 


± 


0.02 


FUCD0336-3514 


19 280 


± 


700 


-0.04 


± 


0.02 


> 15 000 


-0.23 




0.01 


FUCD0336-3548 


7 670 


± 


250 


-0.55 


± 


0.02 


7 980 ± 210 


-0.70 




0.02 


UCDl 


9 710 


± 


260 


-0.06 


± 


0.01 


13 760 ± 350 


-0.42 




0.01 


FUCD0337-3536 


5 400 


± 


170 


-0.60 


± 


0.02 


7 350 ± 150 


-0.87 


± 


0.02 


FUCD0337-3515 


8 410 


± 


260 


-0.06 


± 


0.02 


10 750 ± 250 


-0.49 


± 


0.02 


FUCD0338-3513 


9 530 


± 


580 


-1.20 


± 


0.03 


6 236 ± 99 


-1.16 


± 


0.02 


FUCD0339-3519 


13 000 


± 


1 100 


-0.32 


± 


0.02 


12 980 ± 340 


-0.46 


± 


0.01 


UCDS 


12 450 


± 


280 


-0.97 


± 


0.03 


7 210 ± 160 


-1.01 


± 


0.02 


NGC1407GC1 


6 760 


± 


340 


-0.47 


± 


0.03 


6 290 ± 200 


-0.67 


± 


0.02 


VUCD1230+1233 


11380 


± 


830 


-0.33 


± 


0.02 


9 730 ± 230 


-0.76 


± 


0.02 


VUCD2 


5 473 


± 


87 


-0.49 


± 


0.01 


9 950 ± 200 


-0.98 


± 


0.01 


VUCD3 


24 600 


± 


2 300 


-0.28 


± 


0.09 


> 15 000 


-0.15 


± 


0.01 


Strom417 


4 090 


± 


170 


-0.05 


± 


0.01 


8 500 ± 150 


-0.61 




0.01 


VUCD1231+1234 


14 080 


± 


110 


-1.33 


± 


0.03 


> 15 000 


-1.67 


± 


0.02 


VUCD5 


14 140 


± 


690 


-0.30 


± 


0.01 


12 240 ± 220 


-0.57 


± 


0.01 


VUCD8 


12610 


± 


630 


-1.03 


± 


0.02 


8 550 ± 160 


-1.08 


± 


0.01 


VUCD9 


4 520 


± 


130 


-0.41 


± 


0.01 


6 812 ± 98 


-0.87 


± 


0.01 


VUCD1232+0944 


55 300 


± 


8 800 


-0.83 


± 


0.05 


12 830 ± 800 


-1.02 


± 


0.02 


VUCD1233+0952 


6 260 


± 


270 


-0.37 


± 


0.02 


10 910 ± 450 


-0.99 


± 


0.02 



" Age-metallicity confidence maps for these results can be found in Appendix [XI 
^ The spectral fitting reached the upper age limit of the model grid. 



FUCD0336-3514=ucd3 29.7 (IChilingarian et al.l |2008| ) and 
FUCD0339-3519=F-23 (|Chilingarian et al.l.l201ll l. For three 
of them our age and metallicity estimates agree very well 
with the published data but have much better precision. 
For UCD5, our present estimate falls within the published 
2-(T uncertainty level which can be explained by the low 
metallicity of this object making it a subject to a strong 
age-metallicity degeneracy and hence hampering the deter- 
mination of its stellar population properties in the short 
wavelength range of FL AMES /Giraffe data analysed by 
IChiUngarian et all l|2008h . 



4 DISCUSSION 

In this section we discuss the implications of the age and 
metallicity measurement presented above in terms of forma- 
tion theories for UCDs. Both the Lick index analysis and 
the direct model fitting reveal the UCDs to be an old pop- 
ulation with super-Solar abundances and a wide range of 
metallicity. More generally, we also compare the properties 
of the UCDs to a wide range of other stellar systems. Almost 
without exception they have very high metallicities for their 
luminosities. 

Starting with the Lick index analysis, we found that 
the UCDs in both clusters had generally high (super-Solar) 
alpha element abundances, implying sh ort formation times 
l|Matteuccil . Il994l : IXhomas et al.l . 120051 ). Similar values are 
observed in most globular cluster populations, as distinct 
from present-day galaxies (both dwarf and giant) which have 
lower abundances, indicative of some continuing star forma- 
tion. At first sight, this observation could be taken as evi- 
dence against the threshing model in that the UCDs have 



higher abun dances than dwarf e lliptic al galaxies. However, 
as noted by lEvstigneeva et al.l (|2008l l. if the UCDs were 
stripped from dwarf galaxies at an early epoch, and their 
gas was removed at the same time, this would prevent fur- 
ther star formation, also leading to high abundances. 



The generally old ages of the UCDs also distin- 
guish them from present-day (nucleated) dwarf elliptical 
galaxies. This too argues against recent disruption, but 
not against an early disruption. The spread of age and 
metallicity of the UCDs is consistent with that observed 
for globular clusters, notably those in the Virgo Cluster. 
These results are generally consistent with earlier work, 
but our increased sample has revealed a larger range 
of metallicity. We find metallicities as low as —1.6 dex 
(one in each cluster: VUCD1231-H234 and FUCD0336- 
3536) according to the model fit s . Thi s is in contrast to 
IChilingarian. Cavatte fc Bergondl (|2008t ) who measured a 
smaller UCD sample, finding them all to be metal-rich 
([Fe/H]> —1), although the larger sample measured by 
IChilingarian et al] (|201ll ) included three objects with lower 
metalicity (-1.39 <[Fe/H]< -1.2). We are now finding that 
UCDs exhibit the full range of (low) metalicities shown by 
globular cluster populations. 



Our data also allow us to compare the UCD popula- 
tions between the Virgo and Fornax clusters as we used 
very similar Gemini (North and South) data for both clus- 
ters and identical analysis. Furthermore the objects observed 
had the same luminosity distribution in both clusters. Us- 
ing the model fits we obtain mean ages of (10.4 ± 1.1) Gyr 
for Fornax and (11.2 ± 0.9) Gyr for Virgo with no sig- 
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nificant differences in the distribution^ Similarly there is 
no significant difference in the mean metalicities of [Fe/H] 
= -0.79 ± 0.14 for Fornax and [Fe/H] = -0.90 ± 0.14 
for Virgo. (In each case, for age and metallicity we quote 
the standard error of the respective me an value.) Th i s con- 
tradicts the earlier suggestions made ( Mieske et al.l . l2006l : 
iFirth. Evstigneeva fc Drinkwateij . |200£ ) that the Fornax 
UCD population was younger than that of Virgo due to 
differences in the cluster environment. The cluster environ- 
ments are different (Virgo is larger and much less relaxed), 
but UCDs are mostly only found close to the central galax- 
ies of each cluster, so we would argue that they experience 
similar local environments in the two clusters, resulting in 
similar internal properties. 

For two of the Virg o cluster objects w i th HS T imaging, 
VUCD 3 and VUCD 5, lEvstigneeva et all (|2007l ) published 
dynamical masses and V-band dynamical mass-to-light ra- 
tios: (M/L)dyn = 5.4±0.9 and 3.9±0.6 in Solar units. Using 
our precise age and metallicity measurements, we can de- 
rive their stellar masses and hence estimate their dark mat- 
ter fractions defined as ((M/L)dyn - (M/L),)/(M/L)dyn. 
The y-band stellar mass-to-light ratios for VUCD 3 and 
VUCD 5 are (M/L). = 4.5 ± 0.2 (8.1 ± 0.4) and 2.85 ± 0.13 
(5.0±0.2) [M / L)q^v respectively, wh ere values in parenthe- 
ses correspond to the lSalpeted (| 19551 ) IMF. The correspond- 
ing dark mass fractions are: 15 ± 25 ( — 50 ± 25) and 25 ± 20 
(—28 ±20) per cent. Negative values indicate that the stellar 
mass estimate exceeds the dynamical one, suggesting that 
the Salpeter IMF is not compatible with the observations 
(see discussion bv lChilingarian et al.l . [2"oilh . 

The fact that we do not detect dark matter in these two 
Virgo UCDs is in contrast to previous work suggesting that 
the Virgo UC Ds had higher mass-t o-light ratios than the 
Fornax UCDs. [Ha segan et al. 1 (|2005h reported mass-to-light 
ratios between 6 and 9 for their "probable" UCDs, signifi- 
cantly higher than comparison measurements for UCDs in 
the Fornax Cluster (1-3; see their Figure 11). Although we 
do not have data for a large sample of Virgo UCDs, VUCD 3 
and VUCD 5 are the most massive UCDs in Virgo and for 
these at least we find no evidence of dark matter. 

In Fig. [S] we plot the model-fitted UCD metallicities 
against their luminosities compared to previous measure- 
ments of UCDs and dwarf galaxies, as well as reference 
data for globular clusters and elliptical galaxies. (We use 
the model metallicities rather than the those derived from 
the Lick indices as they use more of the spectra dalta and 
have smaller statistical uncertainties.) It is immediately ob- 
vious that our data are much more precise than most of the 
previous work. With this precision, the two UCDs with lower 
metallicity stand out as being significantly different to the 
rest of the UCD population. We should note, however, that 
the nature of our sample selection means that we cannot ex- 
clude future measurements of intermediate objects. In fact, 
with increasing measurements we now see there is a con- 
tinuous range of objects with intermediate properties (both 
luminosity or metallicity) filling the space between UCDs 
and compact elliptical ("cE" or "M32-like") galaxies. 

The UCDs in Fig. [5] with new Gemini measurements 



* The Kolmogorov-Smirnov, t-tcst and F-tcst all gave p > 0.5 
that the samples were from the same distributions. 



show no evidence of a metallicity-luminosity relation. We 
also measured the metallicity of the objects as a function 
of projected distance from the nearest large galaxy in both 
clusters and found no correlation in those parameters. In 
particular the two objects from the intra-cluster Virgo field 
(VUCD 1232-^0944 and VUCD 1233-^0952 at 550 kpc from 
M87) have metallicities entirely consistent with the other 
Virgo objects. 

Fig. [5] does, however, show that (with two exceptions) 
the UCDs lie well above the canonical metallicity-luminosity 
trend for early-type galaxies (including nucleated dwarf el- 
liptical galaxies) , as was demonstrated previo usly for For- 
nax Cluster UCDs (|Chilingarian et al.l . l20Tll ). This is ex- 
actly what we might expect from tidal disruption of current- 
day dwarf galaxies: the luminosity would decrease while the 
metallicity would remain high. On the other hand, the fig- 
ure also shows that many globular clusters (albeit at lower 
metalicities and luminosities) also lie above the canonical 
metallicity-luminosity trend. Thus the UCDs could repre- 
sent either the high-luminosity end of the globular clus- 
ter metallicity-luminosity relation, or the natural end point 
of a constant-metallicity stripping process of dwarf galax- 
ies. There are now sufficient intermediate objects on both 
sides that the data shown in this figure cannot rule out ei- 
ther hypothesis. A promising approach to help separate the 
two populations would be to measure the ph ysical sizes of 
the ob jects. This has recently been done by iBrodie et al.l 
(|201ll ) who found evidence that UCDs in the Virgo Cluster 
formed a distinct population with significantly larger sizes 
than globular clusters. 



5 SUMMARY 

For each of our Fornax cluster and Virgo cluster UCDs, 
we measured both the Lick/IDS line strength indices and 
performed the nbursts full spectral fitting technique, and 
compared these results to SSP models. Our findings are sum- 
marised below: 

• While much of the evidence, such as super-Solar a- 
abundance ratios, old age measurements, and a large spread 
in the metallicity measurements, seems to imply that the for- 
mation of UCDs is consistent with that of GCs, it does not 
rule out the threshing formation model. The data are only 
inconsistent with stripping of present-day nucleated dwarf 
galaxies, but not with much older disruptions before the nu- 
clei evolved to their current composition. 

• Many of our measurements contrast those from previous 
work. The differences are outlined below: 

(i) Our metallicity measurements exhibit a larger range 
of metallicity than indicated in previous work, consistent 
with the full range of metallicities of GCs. 

(ii) Our age and metallicity measurements, along with lu- 
minosity distributions, indicate similar UCD populations 
in both the Virgo and Fornax clusters. Previous work sug- 
gests that the different cluster environments cause the 
Fornax UCD population to be younger than that of Virgo. 

(iii) We do not detect dark matter in the two Virgo UCDs 
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Figure 5. The metallicity-luminosity relation for UCDs and other stellar systems. The model-fitted metalicities of each UCD are plotted 
a gainst their absolute B ma gnitudes. The dashed line traces the canonical metallicity-luminosity relation for early-type galaxies (e.g. Fig. 8 
of lChilingarian et al.l, I2OIII ') . Comparison values f or diff erent types of galaxy and globular cluste r s are taken fr om the literature accordin g 
to the f ollowing code s : A49 6 - IChihngarian et al.l 1I2OO8I) COS - IChilinga.rian. Cavatte fc Bergondl ||2008|V C09c -IChilingarian et all l[20oi), 
C09v - IChilingarianI l|2009f l. Cll - IChilingarian et a'iri|201ll ), ^98 -iMateJl 

iPrice et al. I 120091). Data for the nearby cE galaxies ar e from IChilinga rian fc Mamoiil 1 120081 ) (M59cO), IChilingarian fc Bergondl 1I2OIOI ') 
('N5846cE'). ISanchez-Blazauez et all l|2006bl 'l CN4486B'I. ICraham (2002 ) (M32 luminosity) and lWorthevI l|2004h fM32 metallicitv) . 



with published dynamical masses and dynamical mass-to- 
light ratios, which happen to be the most massive of the 
Virgo UCDs measured, despite previous work indicating 
that the Virgo UCDs had higher mass-to-Ught ratios than 
the Fornax UCDs. 

• The UCDs did not conform to a metallicity-luminosity 
relation, but did mostly lie above the metallicity-luminosity 
trend for early-type galaxies, which would be consistent with 
formation from present-day tidal stripping of dE,Ns. Inter- 
mediate objects on either side of the data indicate that nei- 
ther hypothesis can be ruled out in this way. 

• The current data provide evidence for both formation 
hypotheses, and are insufficient to rule out either one. 
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APPENDIX A: FIGURES FROM THE 
SPECTRAL FITTING PROCESS 
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Figure A2. NBURSTS spectra, their best-fitting pegase.hr templates jLe Borgne et all . |2004) . fitting residuals and confidence levels of 
the age and metallicity determinations (inner panels) for the Fornax UCDs. 
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Figure A3. NBURSTS spectra, their best-fitting pegase.hr templates jLe Borgne et al] . |2004) . fitting residuals and confidence levels of 
the age and metallicity determinations (inner panels) for the Virgo UCDs. 
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Figure Al. NBURSTS spectra , their best-fitting pegase.hr tem- 
plates l|Le Borgne et aU . l2004h . fitting residuals and confidence 
levels of the age and metallicity determinations (inner panel) for 
the object NGC1407GC1. 
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